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ABSTRACT 


This  research  focused  on  the  experimentation  and  modeling  of  a  series  of  micro-  and 
nanoparticles  in  a  variety  of  media  and  the  response  of  those  materials  to  terahertz  (THz) 
electromagnetic  waves. 

A  series  of  experiments  tested  the  response  of  aluminum  microspheres,  nickel 
microspheres  and  carbon  nano  tubes  (CNTs)  in  epoxy,  liquid  tape  or  SPR955  positive 
photoresist  to  THz  radiation  from  5  THz  to  18  THz.  The  experimental  results  showed  the 
absorption  characteristics  of  the  micro-  and  nanoparticles,  as  well  as  the  media  used.  The 
Transfer  Matrix  Method  was  used  to  obtain  values  for  the  real  and  imaginary  portions  of 
the  refractive  index  of  the  photoresist  to  be  used  in  the  finite  element  modeling  of  the 
composites. 

A  series  of  finite  element  models  was  built  using  the  experimentally  derived  data 
and  simulated  from  5  THz  to  18  THz.  The  models  tested  various  geometries  and  numbers 
of  particles  as  well  as  different  thicknesses  of  the  media.  The  modeling  determined  that 
the  absorption  characteristics  seen  experimentally  were  a  function  of  geometry  and 
particle  size.  Accurate  models  can  be  built  of  aluminum  micro  particles  in  photoresist  for 
the  frequency  spectrum  between  10  and  18  THZ. 
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I.  INTRODUCTION 


A.  OVERVIEW 

Particles  in  the  nano-  and  micro-scale  are  of  great  interest  to  researchers  in 
academia  and  in  the  commercial  realm.  As  technology  continues  to  evolve,  it  becomes 
easier  to  characterize  the  mechanical  properties  of  the  particles  and  their  interaction  with 
other  materials.  The  Navy  is  one  entity  that  is  also  focusing  a  large  amount  of  research 
effort  into  incorporating  nanoparticles  such  as  carbon  nano  tubes  (CNTs)  into  many 
different  applications  such  as  sensors,  armor,  composites,  and  explosives. 

Another  realm  of  technology  that  is  receiving  a  renewed  interest  is  the  part  of  the 
electromagnetic  spectrum  between  microwaves  and  the  infrared  region,  known  as  the  far 
infrared  or  the  terahertz  (THz)  range,  which  is  between  100  GHz  and  30  THz,  and  seen  in 
Figure  1.  It  has  been  studied  in  fields  such  as  astronomy  and  analytical  science  for  many 
years  in  the  past,  but  due  to  recent  advancements  in  nanotechnology  and  photonics,  THz 
technology  is  seeing  use  in  many  different  applications.  Some  of  these  applications 
include  information  and  communications  technology,  non-destructive  evaluation,  global 
environmental  monitoring  and  ultrafast  computing  [1]. 
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Figure  1  Chart  of  the  Complete  Electromagnetic  Spectrum,  from  [2] 


Since  both  of  these  technologies  are  being  explored  and  used  at  an  increased  rate, 
it  is  important  to  understand  if  there  are  any  important  interactions  that  occur  between 
them  as  it  is  very  likely  that  some  form  of  THz  sensing  will  be  used  on  materials 
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containing  micro-  and  nanoparticles.  These  interaetions  ean  help  identify  any 
vulnerabilities  that  may  present  themselves  as  a  result  of  the  interaetions. 

1.  Carbon  Nanotubes,  Aluminum  Microparticles,  and  Nickel 
Microparticles 

One  of  the  great  diseoveries  in  the  realm  of  nanoteehnology  was  the  diseovery  of 
CNTs.  This  diseovery  is  largely  eredited  to  Sumio  lijima  in  1991  in  a  letter  written  to  the 
magazine  Nature  [3].  Sinee  then,  mueh  researeh  has  been  eondueted  in  the  field  of  CNTs. 
Current  researeh  at  the  Naval  Postgraduate  Sehool  (NPS)  is  foeused  on  using  various 
nanopartieles  and  nanotubes  for  applieations  in  energy  generation,  eomposites,  sensors 
and  armor. 

Some  speeifie  examples  are  the  use  of  CNTs  for  the  eonstruetion  and  the  repair  of 
ships  as  part  of  eomposite  patehes  [4],  [5].  Composites  that  are  epoxy-based  are 
frequently  utilized  in  many  applieations  due  to  their  qualities  sueh  as  high  modulus,  high 
ereep  resistanee,  low  shrinkage  at  elevated  temperatures,  high  speeifie  strength  and 
superior  eorrosion  resistanee  [6].  The  researeh  shows  that  ineorporating  CNTs  inereases 
the  values  of  some  of  these  properties.  Other  researeh  has  foeused  on  the  use  of  CNTs  in 
epoxy-based  adhesives  to  hold  together  eomposite  struetures  [7]. 

Other  physieal  properties  of  CNTs  have  led  to  their  use  in  sensors.  They  have 
large  surfaee  area  [8]  and  ean  earry  the  highest  eurrent  density  of  any  metal,  10^-10^° 
A/em^  [9].  This  is  important  when  using  them  in  an  eleetrieal  eontext.  Current  researeh  at 
NPS  is  investigating  the  use  of  earbon  nanofibers  (CNFs)  for  their  use  in  supereapaeitor 
applieations  [10]. 

Aluminum  nano-  and  miero-partieles  have  also  been  used  in  eomposite  struetures. 
Depending  on  the  size  of  partiele  used,  different  meehanieal  properties  are  affeeted.  If 
nano-sized  partieles  are  used,  then  there  is  an  inerease  in  the  work  of  fraeture.  If  mieron- 
sized  partieles  are  used,  there  is  a  greater  inerease  in  the  fraeture  toughness  of  the 
material  [6].  Other  studies  have  shown  that  the  addition  of  mieron-sized  niekel  partieles 
to  a  mixture  of  epoxy  with  aluminum  miero-partieles  further  inereases  properties  sueh  as 
Young’s  Modulus  and  Yield  Strength  [11]. 
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2.  Terahertz  Technology 

Many  materials,  such  as  plastic,  paper,  and  cloth,  are  opaque  in  the  infrared  and 
visible  regions  of  the  electromagnetic  spectrum  are  transparent  to  THz  radiation  [12],  The 
THz  range  has  advantages  versus  X-rays,  when  used  in  sensing  applications,  due  to  the 
fact  that  THz  waves  are  non-ionizing  and  harmless  to  human  tissue.  This  had  led  to  the 
increase  in  the  use  of  THz  in  various  scanning  applications. 

The  use  of  the  THz  range  of  the  electromagnetic  spectrum  has  required  extensive 
research  into  the  materials  that  are  suitable  for  cameras  and  sensors.  Nano-scale  metal 
films  have  studied  extensively  to  determine  the  characteristics  required  for  high 
absorption  in  the  THz  region  [13],  [14].  Sensors  utilizing  these  technologies  can  then  be 
built  and  used  in  applications  involving  security,  biological,  medical,  and  pharmaceutical 
sciences. 

B,  THESIS  OBJECTIVES 

The  goal  of  this  thesis  research  is  to  build  upon  the  knowledge  of  micro-  and 
nanoparticles  and  THz  electromagnetic  waves,  specifically  the  interaction  between  the 
two.  As  the  use  of  THz  electromagnetic  waves  for  sensors  increases,  it  is  important  to 
know  what  effects  it  will  have  on  the  materials  that  are  subjected  to  this  radiation. 
Materials  such  as  composites  that  contain  nanoparticles  show  that  they  provide  favorable 
increases  in  strength  and  toughness,  but  what  is  the  effect  to  absorption  of  THz  waves? 
This  can  lead  to  unknown  heating  effects  within  the  material,  which  over  time  can  have 
an  effect  on  the  properties  of  the  material. 

By  creating  a  series  of  experiments  to  test  the  response  of  various  micro-  and 
nanoparticles  in  a  variety  of  media  to  THz  electromagnetic  waves,  the  reflection, 
transmission,  and  absorption  characteristics  will  be  better  understood.  The  data  obtained 
from  those  experiments  will  be  used  in  the  development  of  finite  element  models  that 
accurately  represent  the  experimental  results.  Those  models  can  then  be  analyzed  to 
understand  the  heating  effects  of  the  radiation  on  both  the  micro-  and  nanoparticles  and 
the  media  that  contain  them. 
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II.  EXPERIMENTAL  SAMPLE  PREPARATION 


The  goal  of  this  research  was  to  determine  the  characteristics  of  various 
nanoparticles  and  their  response  to  THz  electromagnetic  waves  and  to  build  an  accurate 
finite  element  model  of  the  composites  filled  with  micro-  and  nanoparticles.  In  order  to 
accomplish  this,  data  was  collected  from  a  series  of  prepared  samples.  Many  different 
materials  were  mixed  with  a  variety  of  media  and  then  were  subjected  to  THz 
electromagnetic  waves.  This  chapter  details  the  preparation  of  all  of  the  various  samples 
and  the  molds  used  to  make  the  samples. 

A,  SAMPLE  PREPARATION:  MOLDS  USED 

A  variety  of  media  were  used  throughout  the  conduct  of  this  research  to  hold  the 
variety  of  particles  tested.  The  media  used  were  a  two-part  epoxy,  Gardner  Bender  liquid 
tape  and  SPR955  positive  photoresist.  These  media  were  selected  based  on  their 
availability  in  the  laboratory  setting  at  NPS,  their  ease  of  use,  and  there  estimated 
response  in  the  THz  range.  When  the  samples  were  prepared,  a  “blank”  sample  of  each 
medium  was  made  in  order  to  determine  the  THz  absorption  characteristics  of  the 
medium  without  any  particles  in  the  sample. 

Due  to  the  varying  characteristics  for  each  medium,  such  as  viscosity  and  drying 
time,  different  methods  of  sample  preparation  were  used  throughout  the  duration  of  the 
research.  An  aluminum  mold  was  used  for  the  epoxy  samples,  and  silicon  wafers  in 
various  states  were  used  with  the  liquid  tape  and  the  photoresist  samples. 

1,  Epoxy-Based  Substrate  Mold 

An  aluminum  mold,  seen  in  Figure  2,  was  used  to  make  the  epoxy  based  samples. 
The  mold  was  15.24  cm  (6  in)  in  length,  5.08  cm  (2  in)  in  width,  and  1.27  cm  (0.5  in) 
thick.  There  were  two  holes,  5.08  cm  (2  in)  in  diameter  and  2.0  mm  deep  drilled  into  the 
mold  to  form  the  samples.  The  size  of  two  inches  diameter  was  selected  based  on  the  size 
of  the  sample  tray  of  the  MappIR  Accessory  of  the  Fourier  Transform  Infrared 
Spectroscopy  (FTIR),  which  will  be  discussed  later  in  this  thesis.  A  small  hole  was 
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drilled  in  the  center  of  each  mold  to  allow  high  pressure  air  behind  the  sample  in  order  to 
help  remove  it  from  the  mold  when  it  had  dried.  The  NPS  machine  shop  built  the  mold 
for  the  research. 


2"  diameter,  2  mm  deep 


Figure  2  Aluminum  Mold  Used  For  Epoxy  Based  Samples 


Various  methods  were  attempted  in  order  to  determine  the  best  one  to  remove  the 
epoxy-based  samples.  Three  different  sample  release  conditions  of  the  mold  were  used 
when  attempting  to  remove  the  samples.  The  first  method  consisted  of  using  multiple 
coats  of  an  epoxy-release  wax  within  the  holes  of  the  mold.  The  wax  used  was  Meguiar’s 
Mirror  Glaze  88  Universal  Mold  Release  Wax.  The  second  method  consisted  of  using 
aircraft  flash  tape.  This  tape  is  used  when  making  composite  parts  for  aircraft  to  form 
clean  edges  on  the  part  being  made.  The  third  method  consisted  of  using  the  flash  tape 
coated  with  the  release  wax. 

In  addition  to  the  mold  release  conditions,  a  combination  of  high  pressure  air  and 
heat  were  used  when  trying  to  remove  all  of  the  samples.  Initially  the  mold  was  heated 
with  a  propane  torch  in  order  to  release  the  sample.  This  method  was  not  ultimately  used 
for  all  of  the  samples  analyzed  in  the  FTIR  as  the  heat  may  have  changed  some  of  the 
physical  properties  of  the  epoxy-based  samples.  The  best  method  to  remove  the  samples 
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was  to  coat  the  mold  with  at  least  six  eoats  of  release  wax,  allow  the  sample  to  dry 
overnight,  and  the  use  high  pressure  air  to  remove  the  sample. 

2,  Silicon  Wafer  Mold  Used  with  Photoresist  Solution 

A  mold  was  eonstructed  for  media  that  start  out  as  a  fluid  and  then  dry  into  a 
solid.  In  our  case,  this  fluid  was  SPR955  positive  photoresist.  A  300  gm  thiek  silieon 
wafer  was  used  a  base.  A  separate  silieon  wafer  was  cleaved  into  smaller  pieees  sueh  that 
two  windows  were  formed.  These  windows  served  as  the  areas  to  put  the  materials  to  be 
tested.  The  windows  were  approximately  2.4  em  by  3.4  cm,  and  2.1  cm  by  3.5  cm.  The 
exact  size  of  the  windows  was  not  a  major  coneern.  The  windows  needed  to  be 
sufficiently  large  to  have  a  sample  that  could  easily  be  analyzed  in  the  FTIR.  The  cleaved 
pieces  were  attaehed  to  the  silicon  wafer  base  using  Scotch  tape.  The  mold,  with  SPR 
samples  in  it,  is  in  Figure  3.  The  separate  windows  keep  the  photoresist  in  a  single  area 
while  it  dries.  Unlike  the  aluminum  epoxy  mold,  this  mold  is  a  one-time  use  mold. 
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3. 


Silicon  Wafer  Sandwich  Mold 


A  reusable  mold  made  of  three  silieon  wafers  was  construeted  to  hold  various 
types  and  sizes  of  nanopartieles  that  were  not  submerged  in  media  that  were  liquid  or 
solid,  but  in  air.  It  was  made  sueh  that  it  eould  be  fdled  and  emptied  multiple  times,  as 
long  as  the  silieon  wafers  were  not  broken. 

The  mold  was  made  of  one  silieon  wafer  as  a  base.  The  next  layer  was  made  of  a 
silieon  wafer  that  was  broken  into  multiple  seetions  as  to  form  a  window,  where  any 
number  of  nanopartieles  ean  be  held.  Two  windows  were  built  into  this  middle  layer  in 
order  to  maximize  the  number  of  measurements  that  eould  be  taken.  The  windows  are 
approximately  1.5  mm  by  1.6  mm  and  1.1  mm  by  1.5  mm.  Again,  the  size  of  the  window 
only  had  to  be  large  enough  to  be  analyzed  in  the  FTIR.  The  final  layer  is  made  of  a 
eomplete  silieon  wafer.  All  of  the  silieon  wafers  are  approximately  500  pm  thiek.  This 
mold,  without  the  top  silieon  wafer,  is  in  Figure  4.  All  layers  of  this  mold  are  held 
together  with  Seoteh  tape  along  the  edges  of  the  silieon  wafers. 


Figure  4  Silieon  Sandwieh  Mold  for  Various  Nano-  and  Mieropartieles 

B.  SAMPLE  PREPARATION:  MEDIA 

Multiple  media  were  used  in  the  preparation  of  the  samples  for  the  experimental 


portion  of  this  researeh.  Epoxy,  Gardner  Bender  Liquid  Tape,  SPR955  positive 
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photoresist,  and  air  were  all  used  as  media  to  hold  the  various  miero-  and  nanopartieles, 
and  they  are  deseribed  in  the  following  seetion. 

1.  Epoxy  as  a  Medium 

Epoxy  was  used  as  a  medium  for  earbon  nanotubes,  as  well  as  niekel 
mieropartieles.  The  epoxy  used  was  a  two-part  epoxy  made  by  Pro-Set.  The  eomponents 
of  the  epoxy  were  Ml 002  Resin  and  226  Hardener.  This  partieular  epoxy  was  seleeted  as 
a  eandidate  for  this  researeh  based  on  its  availability  in  the  laboratory  setting  at  NPS,  its 
use  in  other  researeh  projeets  at  NPS,  and  its  ease  of  use. 

The  ratio  of  the  resin  to  the  hardener  was  three  parts  resin  to  one  part  hardener. 
For  the  aluminum  mold  diseussed  previously,  approximately  18  mL  of  resin  mixed  with 
6  mL  of  hardener  was  used  to  fill  both  of  the  holes  in  the  mold.  The  mold  was  filled  so 
that  the  epoxy  and  the  epoxy  with  nanotubes/mieropartieles  eombination  was  level  with 
the  top  of  the  mold.  The  samples  were  left  to  air  dry  overnight.  Multiple  samples  made 
only  of  epoxy  were  produeed  to  use  as  eontrol  samples  when  analyzed  with  the  FTIR. 

Multiwall  CNTs  were  used  in  the  epoxy  substrate.  The  eharaeteristies  of  the 
CNTs  ean  be  seen  in  Table  1.  They  were  mixed  with  the  epoxy  in  eoneentrations  of 
0.1%,  1.0%,  4.0%,  and  10.0%  by  weight,  seen  in  Figure  5.  These  weight  pereentages 
were  seleeted  based  on  diseussions  with  other  thesis  researeh  being  eondueted  at  the  time. 
Additionally,  the  amount  of  CNTs  inereased  in  small  inerements  to  see  what  ehanges  in 
absorption,  if  any,  would  be  seen  by  varying  the  pereentages.  As  the  pereentage  of  CNTs 
inereased  to  10.0%,  the  epoxy  mixture  beeame  more  diffieult  to  mix  and  did  not  dry  as 
smoothly  as  the  samples  that  eontained  less  CNTs. 


Table  1  Multiwall  CNTs  Charaeteristies  Used  in  Epoxy  Samples 


Outer  Diameter 

20-40  nm 

Length 

10-30  pm 

Purity 

>95%  wt 

Ash 

<1.5%  wt 
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Pure  epoxy  0.1%  CNT  1%CNT  4%  CNT 

Figure  5  Epoxy-Based  Samples  with  Various  Percentages  of  CNTs 

Nickel  microparticles  were  also  used  with  epoxy  as  a  medium.  The  nickel 
nanoparticles  were  listed  to  be  a  nominal  size  of  5  |um  in  diameter.  A  Zeiss  scanning 
electron  microscope  (SEM)  image  of  the  nickel  powder  was  taken.  Figure  6  shows  that 
the  particles  were  not  all  5  pm  in  diameter,  and  in  fact  many  of  them  are  much  smaller 
than  5  pm. 


(a)  (b) 

Figure  6  SEM  Images  of  Nickel  Powder  with  Nominal  5  pm  Diameter 

Spheres  at  (a)  500x  Magnification  and  (b)  4000x  Magnification 

The  nickel  particles  were  mixed  with  the  epoxy,  maintaining  the  three  to  one  ratio 
of  resin  to  hardener.  A  4.0%  concentration,  by  weight,  of  nickel  nanoparticles  was  added 
to  the  epoxy  and  the  sample  can  be  seen  in  Figure  7. 
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Figure  7  Epoxy-Based  Samples  with  Nickel  5  pm  Diameter  Nanoparticles 

2,  Liquid  Tape  as  a  Medium  with  a  Silicon  Wafer  as  a  Substrate 

A  substance  known  as  liquid  tape  was  selected  as  a  medium  to  suspend  the 
microparticles.  Liquid  Tape  is  a  product  produced  by  Gardner  Bender.  It  is  an  extremely 
viscous  material  typically  used  to  insulate  and  splice  electrical  wires.  When  the  liquid 
tape  dries,  it  becomes  a  solid.  This  material  was  selected  in  an  effort  to  have  a 
homogenous  mixture  of  nanoparticles  throughout  its  thickness.  The  liquid  tape’s  response 
to  THz  electromagnetic  waves  was  unknown  when  it  was  selected.  It  was  believed  that 
the  metal  nanoparticles  would  not  cluster  in  the  bottom  of  the  mixture  since  the  liquid 
tape  dried  quickly. 

A  control  wafer  using  only  liquid  tape  was  made.  The  liquid  tape  was  applied 
with  a  pipette  to  a  300  pm  thick,  four  inch  silicon  wafer  on  a  spin  coater.  The  liquid  tape 
dried  very  quickly,  but  using  the  spin  coater,  it  covered  a  majority  of  the  surface  of  the 
wafer  as  seen  in  Figure  8.  As  expected,  it  dried  quickly  and  was  selected  as  a  candidate  to 
add  the  aluminum  nanoparticles. 
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Figure  8  Liquid  Tape  Applied  with  a  Spin  Coater  Onto  a  300  pm  Thick,  10.16 

cm  (4  in)  Diameter  Silicon  Wafer 

The  liquid  tape  was  mixed  with  20  pm  diameter  aluminum  nanoparticles,  with  a 
concentration  of  20%  by  weight.  An  SEM  image.  Figure  9,  of  the  aluminum  powder  was 
taken  to  verify  the  diameters  of  the  aluminum  spheres.  The  liquid  tape,  with  the 
aluminum  nanoparticles,  was  then  applied  to  a  300  pm  thick  silicon  wafer  using  a  spin 
coater.  The  silicon  wafer  was  selected  as  a  substrate  so  that  the  liquid  tape  and  aluminum 
mixture  would  have  flat  surface  to  adhere  to. 


(a) 


(b) 


Figure  9  SEM  Images  of  Aluminum  Power  with  Nominal  20  pm  Spheres  at 

(a)  lOOx  and  (b)  500x  Magnification 
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Using  liquid  tape  was  difficult  due  to  its  high  viscosity,  especially  after  contact 
with  air.  It  was  difficult  to  get  a  smooth  surface  using  the  spin  coater  because  the  liquid 
tape  had  a  very  rapid  drying  time,  even  with  the  addition  of  a  solvent.  The  uneven 
spreading  of  the  mixture  can  be  seen  in  Figure  10.  To  further  investigate  the  effects  of  the 
uneven  spreading  of  the  liquid  tape  on  the  silicon  wafer,  an  SEM  image  of  a  piece  of  the 
wafer  was  taken.  Due  to  the  amount  of  nonconductive  material  in  the  liquid  tape,  the 
sample  had  to  be  sputter-coated  with  a  small  layer  of  platinum  in  order  to  obtain  a  useful 
SEM  image  of  the  aluminum  particles  suspended  in  the  liquid  tape,  seen  in  Eigure  11.  It 
can  be  seen  that  there  was  some  clumping  of  the  aluminum  particles,  as  many  of  the 
particles  seen  in  the  SEM  image  are  much  larger  than  20  pm.  The  clumping  of  the 
particles  as  well  as  the  natural  distribution  of  uneven  sizes  of  the  aluminum  powder  can 
have  some  significant  effects  in  the  absorption  of  THz  at  different  frequencies. 


Eigure  1 0  Eiquid  Tape  Mixed  with  20  pm  Aluminum  Nanoparticles,  Applied 
with  a  Spin  Coater,  on  a  300  pm  Silicon  Wafer 


13 


(a)  (b) 

Figure  1 1  SEM  Images  of  Sputter  Coated  Liquid  Tape  Sample  with  Nominal 
20  pm  Diameter  Aluminum  Spheres  in  Liquid  Tape  Showing 
Clumping,  Seen  at  (a)  lOOx  and  (b)  500x 


3.  Photoresist 

SPR955  positive  photoresist  was  selected  as  another  medium  for  the 
microparticles.  The  photoresist  was  selected  as  it  is  fairly  transparent  in  the  THz  range 
when  applied  in  a  thin  film.  The  viscosity  of  the  photoresist  is  such  that  particles  are 
easily  mixed  into  the  solution.  The  20  pm  diameter  aluminum  spheres,  discussed 
previously,  were  mixed  into  the  photoresist.  The  photoresist  was  put  into  the  silicon 
wafer  mold,  previously  discussed,  designed  to  be  used  for  this  purpose.  After  the 
photoresist  was  applied  to  the  mold,  it  was  put  onto  a  100  degree  C  hotplate  to  dry  for 
approximately  10  minutes.  The  mold  was  set  up  such  that  one  window  was  for  a  control 
sample  of  only  photoresist.  The  second  window  was  for  the  photoresist  mixed  with  the 
aluminum  particles.  Unfortunately,  it  was  difficult  to  get  an  accurate  measurement  of  the 
amount  of  aluminum  particles  that  were  added  to  photoresist  solution  due  to  the  viscosity 
of  the  fluid  and  the  added  precaution  of  working  within  a  vapor  hood  when  working  with 
the  fluid. 

4.  Various  Nanoparticles  in  Air 

Various  quantities  of  CNTs,  aluminum  and  nickel  microspheres  were  tested  in  the 
silicon  wafer  sandwich-type  mold.  For  these  experiments,  the  material  was  simply  poured 
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into  the  windows  formed  on  the  silieon  wafers  until  the  window  was  completely  filled. 
The  mold  was  weighed  prior  to  and  after  each  type  of  particle  was  added  to  determine  to 
approximate  weight  of  material  added.  Table  2  lists  the  various  nanoparticles  and 
amounts  used  in  this  research.  Figure  12  shows  the  sandwich  mold  when  it  was  filled 
with  the  aluminum  and  the  nickel  particles. 


Table  2  Mass  and  Types  of  Particles  Used  in  Silicon  Wafer 

Sandwich  Mold 


Material 

Mass  Used  (grams) 

Multi-Walled  CNTs 

0.0394 

5  pm  diameter  Nickel  Spheres 

0.0504 

20  pm  diameter  Aluminum  Spheres 

0.1305 

Figure  12  500  pm  Thick  Silicon  Wafer  Mold  with  20  pm  Aluminum  Powder 

and  5  pm  Nickel  Powder 


15 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


16 


III.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 


All  samples  were  analyzed  using  the  Thermo  Scientific  Nexus  870  FTIR  with  a 
Pike  Technologies  MappIR  Automated  Semi-Conductor  Wafer  Accessory  and  a  standard 
desktop  personal  computer.  All  of  the  experimental  results  were  analyzed  from  5  THz  to 
18  THz.  The  results  were  used  to  build  and  fine  tune  the  COMSOL  models  of  the  various 
samples. 

A,  FTIR  OPERATION 

In  order  to  analyze  the  response  of  the  particles  in  various  media  to  THz 
electromagnetic  waves,  an  FTIR,  seen  in  Figure  13,  was  utilized.  The  device  consists  of  a 
Michelson  interferometer  that  has  a  mirror  attached  to  a  moving  arm.  The  moving  arm 
produces  periodic  oscillations  in  amplitude  which  changes  the  frequency  based  on  the 
wavelength  of  radiation  and  the  velocity  of  the  moving  mirror.  There  is  a  thermal 
detector  within  the  FTIR  that  measures  the  frequency  data  and  performs  a  Fourier 
Transform  on  the  data.  This  process  is  performed  across  the  whole  spectrum  of  concern, 
which  for  the  purpose  of  this  research  consists  of  wavelengths  between  16.7  pm  and 
60  pm. 


Figure  13  Thermo  Scientific  Nexus  870  FTIR  Spectrometer 
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The  FTIR  available  for  use  in  the  NFS  SRL  contains  a  globar  source,  consisting 
of  a  hot  SiC  filament,  and  a  Deutrated  Triglycine  Sulfate  (DTGS)  pyroelectric  detector 
with  a  polyethelene  window.  This  combination  of  source  and  detector  is  recommended 
by  the  manufacturer  for  measurements  in  the  far  infrared  spectrum,  between  700  cm‘^  to 
20  cm‘*  [15],  which  equates  to  0.6  THz  to  approximately  21  THz.  This  setup  is  optimized 
for  the  spectrum,  5  THz  to  1 8  THz,  of  concern  for  this  research. 

The  Nicolet  870  is  fitted  with  a  Pike  Technologies  MappIR  Automated  Semi- 
Conductor  Wafer  Accessory.  The  MappIR  Accessory  allows  the  automated  analysis  of 
multiple  points  on  a  standard  4-inch  silicon  wafer  or  the  2-inch  epoxy  samples  utilized  in 
this  research.  The  MappIR  Accessory  uses  AutoProS  software,  seen  in  Figure  14,  to 
control  the  rotation  and  movement  of  the  sample.  The  required  scanning  pattern  is 
programmed  within  AutoPro5  by  the  user.  The  sample  tray  can  rotate  360  degrees  and 
moves  the  sample  to  the  positions  programmed  by  the  user  using  a  radius  and  an  angle  of 
rotation  from  the  starting  position. 


Figure  14  AutoPro5  Interface 

The  FTIR  uses  the  OMNIC  program  for  the  actual  spectral  analysis  and  data 
collection  of  a  sample.  The  OMNIC  program  is  also  used  to  perform  the  alignment  of  the 
laser  after  switching  modes  between  reflection  or  transmission. 

The  FTIR  has  two  modes  of  operation.  One  measures  the  transmittance  of  a 
sample  and  the  other  measures  the  reflectance  of  the  material.  The  absorption,  which  is 
the  ultimate  concern  throughout  the  course  of  this  research,  is  a  calculated  value  where 
Equation  (1)  is  used. 

%Absorption  -100-  %Reflectance  -  %Transmittance  ( 1 ) 
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Changing  the  mode  of  operation  is  done  on  the  MappIR  Aecessory.  The 
configuration  can  be  seen  in  Figure  15.  Prior  to  each  data  collection  run,  whether  it  is  to 
record  reflectance  or  transmittance,  an  alignment  check  and  background  collection  is 
required.  The  alignment  verifies  the  laser’s  signal  prior  to  operation. 


(a)  (b) 

Figure  15  (a)  FTIR  set  up  to  measure  reflection,  from  [16] 

(b)  FTIR  set  up  to  measure  transmission,  from  [16] 

The  background  check  procedure  varies  based  on  whether  the  FTIR  is  set  up  to 
measure  reflectance  or  transmittance.  The  background  check  is  critical  to  the  proper 
operation  of  the  FTIR  as  it  measures  the  100%  signal  value.  For  a  reflective  background 
check,  a  silicon  wafer  coated  with  a  thin  film  of  gold  is  loaded  into  the  FTIR.  The  gold- 
coated  silicon  wafer  acts  as  a  perfect  mirror  that  has  close  to  100%  reflection  in  the  THz 
range.  For  the  transmittance  background  check,  the  FTIR  is  not  loaded  with  a  sample  or  a 
mirror  and  is  run  empty. 

All  samples  were  analyzed  from  approximately  1.50  THz  to  18  THz.  The  data 
obtained  from  1.50  THz  to  5.00  THz  is  not  reliable  due  to  being  near  the  edge  of  the 
beam  frequency  range  with  a  very  weak  signal,  resulting  in  a  large  amount  of  noise  in 
this  range.  Therefore,  this  range  was  omitted  from  all  results. 
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Samples  that  used  the  silicon  wafers  as  substrates  also  encountered  a  phenomenon 
known  as  multiple-beam  interference  in  a  parallel  plate,  or  the  Fabry-Perot  effect,  seen  in 
Figure  16.  An  incoming  wave  of  amplitude.  Eg,  and  at  an  angle  0  reaches  the  first  surface 
of  the  silicon  wafer.  The  silicon  wafer  has  reflection  and  transmission  coefficients  of  r 
and  t,  respectively,  and  internal  reflection  and  transmission  coefficients  of  r’  and  t’, 
respectively.  A  portion  of  the  wave  reflects  off  of  the  surface  of  the  wafer  back  toward 
the  source  with  an  amplitude  of  rEg.  The  rest  of  the  wave  passes  through  to  the  back  side 
of  the  wafer.  Fiere  a  portion  of  the  wave  reflects  back  into  the  wafer,  with  an  amplitude  of 
trEg,  and  a  portion  leaves  the  wafer  with  an  amplitude  of  ttEg.  The  portion  of  the  wave 
that  reflects  within  the  wafer,  having  an  amplitude  of  tr  Eg,  reaches  the  front  side  of  the 
wafer  and  part  of  it  reflects  back  into  the  wafer  and  part  of  it  passes  back  through  the 
front  of  the  wafer.  This  series  of  waves  leads  to  constructive  and  destructive  interference 
patterns  when  the  waves  are  recorded  after  passing  through  the  wafer. 
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Figure  16  Multiple-Beam  Interference  Pattern  (Fabry-Perot  Effect)  Due  to 
Constructive  and  Destructive  Interference,  after  [17] 


Since  the  Fabry-Perot  effect  makes  it  difficult  to  discern  absorption  of  the 
analyzed  film  from  the  effects  of  the  substrate,  a  program  was  written  and  run  in 
MATLAB  to  determine  the  local  maximum  on  the  experimental  absorption  results.  The 
local  maximum  are  then  plotted  and  displayed  for  all  samples  that  have  silicon  wafers  as 
substrates. 
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B. 


RESULTS 


Various  combinations  of  the  previously  described  media  and  micro-  and 
nanoparticles  were  tested  in  the  FTIR.  Those  results  are  contained  in  the  following 
section. 


1,  Epoxy  as  a  Medium 

The  2"  diameter,  2  mm  thick  epoxy-only  samples  were  analyzed  in  the  FTIR.  The 
reflectance  and  transmittance  were  measured  at  three  points  on  the  sample,  seen  in  Figure 
17.  The  sample  points  were  the  center  of  the  sample  and  12  mm  from  the  center  of  the 
sample  at  the  12  o’clock  and  6  o’clock  positions.  All  three  points  had  similar 
characteristics.  Of  note  is  the  fact  that  all  three  points  absorb  90%  to  100%  of  the 
incoming  THz  radiation  between  3  and  18  TFlz.  These  results  are  in  Figure  18.  The  low 
transmission  by  the  epoxy  is  likely  the  cause  of  some  noise  in  the  signal  which 
contributes  to  the  “unclean”  nature  of  the  absorption  results  for  all  samples  with  epoxy. 
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Figure  18  Experimentally  Obtained  Reflection,  Transmission  and  Absorption 
Values  of  2mm  Thick  Epoxy-Only  Sample 


All  epoxy  samples  with  CNTs  were  analyzed  in  the  FTIR.  Due  to  the  similarities 
in  the  results  only  the  epoxy  sample  with  10%  CNTs  (by  weight)  is  discussed  in  this 
manuscript.  The  sample  was  analyzed  at  two  points,  one  at  the  center  of  the  sample  and 
one  4  mm  from  the  center.  The  absorption  of  this  sample  was  close  to  99%  as  seen  in 
Figure  19. 
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Figure  19  Absorption  of  THz  Radiation  in  Epoxy  with  10%  CNTs 

The  absorption  values  obtained  at  each  frequency  value  were  averaged  for  the 
epoxy-only  sample  and  the  epoxy-based  sample  containing  10%  CNTs.  The  averaged 
values  of  absorption  for  the  epoxy-only  sample  were  then  subtracted  from  the  averaged 
absorption  of  the  epoxy  with  10%  CNTs  sample.  The  results  are  contained  in  Figure  20. 
A  simple  linear,  least  squared  average  was  performed  on  the  results  in  the  graph  and  it 
was  found  that  the  CNTs  absorb  approximately  4%  of  the  THz  radiation  between  5  THz 
and  18  THz.  This  is  a  very  small  percentage  of  the  incoming  radiation.  This  small 
number,  make  it  difficult  to  truly  see  the  actual  absorption  due  to  the  CNTs. 
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Figure  20  Absorption  Due  Only  to  CNTs 


The  sample  made  with  epoxy  and  4%  (by  weight)  nickel  particles  was  also  tested 
in  the  FTIR.  The  results,  in  Figure  21,  were  similar  to  the  epoxy  with  10%  CNTs. 
Approximately  94%  of  the  THz  radiation  was  absorbed  by  the  sample.  Again  this  is  due 
to  the  strong  absorption  characteristics  of  the  epoxy  in  the  THz  range  of  the  spectrum. 
The  absorption  of  the  nickel  particles  alone  was  calculated  by  subtracting  the  absorption 
of  the  epoxy-only  sample  from  the  sample  with  nickel  particles,  and  is  in  Figure  21.  The 
absorption  by  the  nickel  particles  is  less  than  the  absorption  by  the  CNTs. 
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Figure  21  Absorption  Characteristics  of  4%  (by  Weight),  5  pm  Diameter 

Nickel  Spheres  in  Epoxy 
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Figure  22  Absorption  Due  Only  to  5  pm  Nickel  Particles 


Overall,  the  results  from  both  the  CNTs  and  nickel  particles  in  epoxy  show  that 
epoxy  exhibits  strong  absorption  in  the  TFIz  range  of  the  spectrum  when  the  thickness  is 
approximately  2  mm.  Due  to  this  fact,  it  is  not  the  best  medium  to  use  when  testing  the 
absorption  characteristics  of  the  various  nanoparticles. 
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2,  Liquid  Tape  as  a  Medium 

As  with  the  epoxy-based  samples,  all  of  the  samples  using  the  liquid  tape  as  a 
medium  were  analyzed  using  the  FTIR  to  obtain  reflection  and  transmission 
characteristics.  For  the  sample  containing  only  liquid  tape  on  a  silicon  wafer  substrate, 
three  sample  points  were  selected  for  analysis.  One  sample  point  was  in  the  center,  with 
the  other  two  points  located  12  mm  from  the  center,  along  the  vertical  centerline  of  the 
sample.  The  calculated  absorption  values  at  the  three  sample  points  were  averaged  and 
plotted  as  seen  in  Figure  23.  There  are  very  distinct  peaks  as  the  absorption  varied 
dramatically  throughout  the  range  of  interest. 

The  sample  containing  liquid  tape  mixed  with  20%  (by  weight)  aluminum 
particles  exhibited  absorption  characteristics  as  seen  in  Figure  24.  The  calculated 
absorption  data  from  Figure  23  and  Figure  24  were  combined  in  the  way  described  in 
previous  chapter  to  extract  the  effect  due  only  to  the  aluminum  particles.  The  combined 
data  is  captured  in  Figure  25. 
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Figure  23  Absorption  Characteristics  of  a  Thin  Layer  of  Liquid  Tape  on  a  300 

pm  Silicon  Wafer 
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Figure  24  Absorption  Characteristics  of  Liquid  Tape  with  20%  Aluminum 

Spheres  on  a  300  pm  Silicon  Wafer 
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Figure  25  Absorption  of  Only  20  pm  Aluminum  Spheres 

3,  Photoresist  as  a  Medium 

FTIR  measurements  were  taken  on  the  sample  prepared  with  the  SPR955  positive 
photoresist  as  the  medium.  The  reflection  and  transmission  were  measured  at  five 
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different  points  on  the  photoresist  only  window.  The  window  that  eontained  a  mixture  of 
photoresist  and  aluminum  partieles  was  measured  at  three  different  points.  .  The  average 
absorption  in  the  photoresist  only  sample  was  ealculated  and  can  be  seen  in  Figure  26. 
The  averaged  absorption  of  the  three  points  of  the  photoresist  with  aluminum 
nanoparticles  sample  is  contained  in  Figure  27.  The  absorption  due  to  the  aluminum 
particles  alone  was  calculated  from  the  averaged  absorption  of  the  photoresist  alone 
subtracted  from  the  averaged  absorption  of  the  photoresist  with  aluminum  particles  and 
can  be  seen  in  Figure  28. 

Again  it  can  be  seen  that  the  absorption  due  to  the  aluminum  particles  alone  is  a 
very  small  number  between  5  THz  and  18  THz,  and  it  is  on  the  same  order  of  magnitude 
as  was  previously  seen  on  the  samples  of  aluminum  in  epoxy  samples.  However,  the 
aluminum  particles  in  epoxy  maintain  a  fairly  consist  absorption  of  4%  through  the  whole 
range  where  the  aluminum  in  photoresist  tapers  off  toward  zero.  This  can  be  due  to  the 
random  size  distribution  of  particles  in  the  two  samples  or  different  amounts  of  aluminum 
in  each  sample. 

One  additional  issue  encountered  during  the  testing  of  the  photoresist  solution 
was  the  fact  that  in  early  testing  the  sample  did  not  seem  to  be  fully  dry.  As  testing 
continued  it  appeared  to  be  more  dry.  This  can  explain  some  of  the  differences  seen  in  the 
absorption  characteristics  of  the  samples  as  the  amount  of  water  in  the  photoresist 
decreased  due  to  evaporation. 
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Figure  26  Averaged  Absorption  of  Photoresist  Only 
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Figure  27  Averaged  Absorption  of  Photoresist  with  20  |um  Aluminum  Particles 
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Figure  28  Averaged  Absorption  of  20  pm  Aluminum  Partieles  Alone 

4,  Various  Particles  in  Air 

The  silieon  wafer  sandwieh  mold  was  analyzed  using  a  variety  of  partieles.  In 
order  to  see  the  amount  of  absorption  due  to  each  type  of  particle,  results  from  the 
absorption  due  to  the  silicon  wafer  and  air  obtained  from  analyzing  the  blank  window  in 
the  sandwich  mold  were  subtracted  from  the  results  from  the  mold  filled  with  the  various 
particles.  Each  window  of  the  silicon  wafer  was  analyzed  at  three  different  points  in  the 
center  of  each  window.  The  points  were  in  a  horizontal  line,  spaced  2  mm  apart. 

The  first  sample  analyzed  was  the  blank  window,  filled  only  with  air.  The 
absorption  values  obtained  at  all  three  sample  points  were  averaged  and  the  results  are  in 
Figure  29.  The  absorption  is  due  to  two  500  pm  thick  silicon  wafers  and  a  500  pm  thick 
airspace  between  the  wafers.  The  absorption  ranges  between  25%  and  70%  due  to  Fabry- 
Perot  effect. 
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Figure  29  Absorption  Characteristics  of  the  Average  of  Three  Sample  Points 
on  the  Silicon  Sandwich  Mold,  Blank  Window 

The  multi-walled  CNT  sample  was  analyzed  at  three  points  and  the  absorption 
values  of  the  three  points  were  averaged,  seen  in  Figure  30.  Those  values  were  then 
subtracted  from  the  averaged  absorption  values  of  the  blank  window  of  the  silicon 
sandwich  mold  in  order  to  determine  the  absorption  due  to  the  CNTs  alone. 
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Figure  30  Averaged  Absorption  Characteristics  of  Multi- Walled  CNTs  in  Air 

in  the  Silicon  Sandwich  Mold 
31 


As  seen  in  Figure  31,  the  multi-walled  CNTs  absorb  close  to  35%  of  the  incoming 
radiation  at  5  THz  but  then  trail  off  to  approximately  7%  close  to  18  THz.  This  is  more 
absorption  than  was  seen  in  the  10%  CNTs  in  epoxy  sample  between  5  THz  and  18  THz. 
The  general  trend  of  the  CNTs  in  air  is  also  different  than  the  CNTs  in  epoxy.  The  CNTs 
in  epoxy  maintain  a  fairly  constant  level  of  absorption  between  5  THz  and  1 8  THZ, 
where  the  CNTs  in  air  had  much  higher  absorption  low  in  the  frequency  band  and  it 
decreases  as  it  approaches  18  THZ.  The  difference  is  likely  due  to  the  difference  in  the 
density  of  CNTs  in  both  instances.  In  the  epoxy,  the  density  of  the  CNTs  was  0.04157 
g/cm  ,  where  the  density  of  CNTs  in  the  air  was  0.32833  g/cm  .  The  volume  density  of 
CNTs  in  the  sandwich  mold  is  almost  ten  times  greater  than  that  of  the  CNTs  in  epoxy. 
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Figure  3 1  Averaged  Absorption  Due  Only  to  Multiwall  CNTs  with  Size 

Characteristics  as  Found  in  Table  1 


The  aluminum  particles  were  analyzed  at  the  same  points  as  the  multi-walled 
CNT  sample  and  the  results  are  in  Figure  32.  Again,  the  values  of  the  blank  window  were 
subtracted  from  the  averaged  values  obtained  from  the  aluminum  particle  analysis.  These 
results,  as  well  as  the  aluminum  particles  in  the  liquid  tape  and  the  photoresist,  are  in 
Figure  33.  Like  the  previous  results  concerning  the  aluminum  powder,  there  is  decreasing 
trend  in  toward  zero  as  the  THz  wave  goes  from  5  THz  to  18  THz.  When  the  absorption 
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due  to  the  aluminum  particles  in  air  is  compared  to  the  aluminum  in  the  photoresist,  the 
aluminum  in  the  air  absorbs  more  of  the  incoming  THz  wave.  This  result  is  expected,  like 
with  the  CNTs,  because  the  concentration  of  aluminum  is  much  greater  for  the  silicon 
wafer  mold  versus  the  aluminum  in  the  photoresist.  The  aluminum  particles  in  the  liquid 
tape  showed  more  absorption  through  the  entire  frequency  spectrum,  even  though  they 
were  only  mixed  with  a  20%  by  weight  ratio  into  the  liquid  tape.  The  relatively  high 
absorption  rate  between  5  THz  and  8  THz  could  be  due  to  the  clumping  of  the  aluminum 
particles  seen  previously  in  this  thesis  as  well  as  presence  of  particles  larger  than  20  pm. 
The  larger,  clumped,  aluminum  particles  absorb  a  greater  amount  of  the  incoming  THz 
wave  low  in  the  spectrum  of  interest.  Another  factor  in  the  high  experimental  absorption 
rates  of  the  aluminum  particles  in  liquid  tape  is  the  fact  that  there  may  be  some  scattering 
effects  that  are  not  accounted  for  in  the  experimental  measurements.  One  assumption 
when  using  the  FTIR  is  that  the  incoming  wave  is  transmitted,  reflected,  or  absorbed. 
There  is  not  a  way  to  measure  any  waves  that  are  scattered  within  the  FTIR. 
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Figure  32  Absorption  of  20  pm  Aluminum  Spheres  in  Air,  in  Silicon  Wafer 

Sandwich  Mold 
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Figure  33  Absorption  Due  only  to  20  pm  Aluminum  Spheres 


Finally,  the  niekel  particles  were  analyzed  using  the  silicon  sandwich  mold  and 
the  results  are  in  Figure  34.  The  same  type  of  analysis  was  performed  as  discussed 
previously  for  CNTs  and  aluminum  to  determine  the  absorption  due  only  to  the  nickel 
particles.  Figure  35  shows  that  the  nickel  particles  absorb  less  than  10%  of  the  incoming 
THz  wave  for  most  of  the  spectrum  and  it  decreases  toward  zero  as  the  values  move  to 
1 8  THz.  This  result  is  similar  to  the  absorption  characteristics  of  the  nickel  particles  in 
epoxy.  The  magnitude  of  absorption  is  higher  for  the  nickel  particles  in  air,  but  like  the 
aluminum  particles,  this  could  be  due  to  the  increased  volume  of  nickel  particles  in  this 
mold. 

In  general,  both  the  nickel  and  aluminum  particles  absorb  a  very  small  percentage 
of  the  incoming  THz  wave  and  both  types  of  particles  show  decreasing  trends  from  5 
THz  to  18  THz. 
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Figure  34  Average  Absorption  of  5  pm  Nickel  Spheres  in  Air,  in  Silicon  Wafer 

Sandwich  Mold 
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Figure  35  Average  Absorption  Due  Only  to  5  pm  Nickel  Spheres 


C.  DETERMINATION  OF  REFRACTIVE  INDEX 

Many  of  the  materials,  such  as  the  epoxy  and  the  photoresist,  used  as  media 
throughout  this  research  did  not  have  readily  accessible  values  for  the  real  (n)  and 
imaginary  (k)  components  of  refractive  index.  Additionally,  both  the  real  and  imaginary 
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components  of  the  refraetive  index  are  frequeney-dependent  at  those  frequeneies.  The 
refraetive  index  values  are  required  inputs  into  the  COMSOL  models  and  are  needed  for 
an  aeeurate  model  simulation.  The  method  used  to  determine  the  refraetive  indiees  takes 
the  experimental  refleetion  and  transmission  measurements  and  applies  the  Transfer 
Matrix  Method.  The  basie  theory  is  deseribed  by  Pedro tti  and  Pedro tti  in  [17].  The 
speeifie  method  used  in  this  researeh  is  deseribed  by  Martin  in  [18]. 

In  order  to  estimate  the  refraetive  indiees  of  the  materials  used  throughout  this 
researeh  a  Thin  Film  Analysis  Tool  was  used.  The  tool  was  ereated  by  Dr.  Fabio  Alves,  a 
professor  at  NPS,  and  Major  Miehael  Martin  developed  the  tool.  The  tool  was  designed 
sueh  that  a  graphieal  user  interfaee  (GUI),  designed  in  Labview  2012,  allows  a  user  to 
adjust  various  parameters  so  that  the  refleetion,  transmission,  and  absorption  eurves 
mateh  the  experimentally  derived  values.  The  Transfer  Matrix  Method  ealeulations  are 
performed  by  MATLAB,  so  the  user  must  have  MATLAB  and  Labview  2012  installed 
on  their  eomputer  for  the  program  to  run  properly. 

The  Fabry-Perot  effeets  did  make  it  more  diffieult  to  preeisely  mateh  the 
experimental  and  simulated  data.  Figure  36  shows  an  example  of  the  simulated  and 
experimental  refleetion,  transmission,  and  absorption  eurves  that  were  required  to  be 
matehed  to  determine  n  and  k.  The  Fabry-Perot  effeets  are  elearly  visible  in  Figure  36. 
After  the  simulated  and  experimental  values  were  matehed,  a  plot  of  n  and  k  versus 
frequeney  was  ereated  to  show  the  frequeney  dependenee  of  those  values,  as  seen  in 
Figure  37  and  Figure  38. 
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Figure  36  (a)  Simulated  Reflection,  Transmission  and  Absorption  in  Transfer 

Matrix  Method  GUI 

(b)  Experimental  Reflection,  Transmission,  and  Absorption  in  Transfer 

Matrix  Method  GUI 
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Figure  37  Thin  Film  Analysis  Values  and  Linear  Fit  for  SPR955  Real  Portion 

of  Refractive  Index 
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Figure  38  Thin  Film  Analysis  Values  and  Second  Order  Polynomial  Best  Fit 
for  SPR955  Imaginary  Portion  of  Refractive  Index 
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IV.  MODELING 


The  software  used  for  the  modeling  portion  of  this  researeh  was  the  COMSOL 
Multiphysies  program,  using  the  RF  module  to  simulate  electromagnetie  waves.  The  RF 
module  solves  Maxwell’s  equations.  COMSOL  was  seleeted  due  to  its  ability  to  model 
multiple  physies  interaetions  within  a  model.  COMSOL  has  been  used  to  support  various 
researeh  topies  involving  the  THz  portion  of  the  eleetromagnetic  speetrum  at  NFS. 


A,  COMSOL  RF  MODULE 


COMSOL  solves  the  series  of  equations  known  as  Maxwell’s  equations.  The 
equations  handle  some  basie  eleetromagnetic  quantities.  The  quantities  are  the  electric 
field  intensity  (E),  the  electric  flux  density  (D),  the  magnetic  field  density  (H),  the 
magnetic  flux  density  (B),  the  current  density  (J)  and  the  electric  charge  density  (p).  The 
form  of  Maxwell’s  equations  used  by  COMSOL  can  be  seen  in  Equations  (2)-(5): 


^  5D 

VxH  =  J+ — 

(2) 

dt 


VxE= - 

(3) 

dt 

V»D  =  p 

(4) 

V.B  =  0 

(5) 

Within  each  COMSOL  model,  the  equations  are  modified  in  order  to  account  for 
conductive  or  dielectric  material  and  they  rely  on  a  series  of  constitutive  relationships. 
The  materials  within  the  model  that  are  conductive  use  Equation  (6)  as  the  governing 
equation  and  Equation  (7)  as  a  constitutive  relationship: 


V  X  (V  X  E)  -  kl  {s^  -  -^)E  =  0 

COSq 

B  = 


(6) 

(V) 
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where  )j,r  is  the  relative  permeability  of  the  material,  Sr  is  the  relative  permittivity  of  the 
material,  a  is  the  conduetivity  of  the  material.  The  required  inputs  to  COMSOL  for 
eonductive  materials  are  pr,  £r,  and  a. 

Materials  that  are  a  dieleetrie  use  Equations  (8)-(ll)  for  their  governing  and 
constitutive  equations: 


Vx(VxE)-ko£^E  =  0 

(8) 

(9) 

cr  =  0 

(10) 

(11) 

where  n  is  the  real  component  of  the  refractive  index  and  k  is  the  imaginary  component 
of  the  refractive  index.  COMSOL  requires  n  and  k  as  material  properties  for  the  materials 
that  are  dielectrics. 

B,  MODELING  A  SINGLE  UNIT  CELL 

One  particular  topic  related  to  this  current  research  was  a  paper  written  by  Alves 
et  al.  As  part  of  the  research  for  that  paper  the  authors  created  a  COMSOL  model  made 
up  of  a  single  unit  cell  of  metamaterial  [19].  This  basic  model  was  modified  in  order  to 
better  fit  this  research.  A  representative  model  used  in  this  research  can  be  seen  in  Ligure 
39.  The  incoming  THz  radiation  is  denoted  by  k,  with  the  associated  magnetic  and 
electric  fields,  shown  by  B  and  E,  respectively. 
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Figure  39  Schematic  of  3D  Finite  Element  Model  of  a  Unit  Cell  Consisting  of 

a  Single  Aluminum  Sphere  in  Liquid  Tape 


The  boundary  conditions  perpendicular  to  the  magnetic  field  are  Perfect  Magnetic 
Conductors.  The  equation  for  this  boundary  condition  is: 

nxH  =  0  (12) 

This  boundary  condition  ensures  that  the  tangential  component  of  the  magnetic  field  is 
zero  and  the  surface  current  is  zero.  This  allows  it  to  be  used  as  a  symmetry  type 
boundary  which  is  useful  in  modeling  a  unit  cell  to  simulate  the  infinitely  repeating  film. 

Similarly,  the  boundary  conditions  perpendicular  to  the  electric  field  are  Perfect 
Electric  Conductors.  This  equation  is: 

nxE  =  0  (13) 

Again,  this  boundary  condition  ensures  that  the  tangential  component  of  the  electric  field 
is  zero  and  is  needed  to  replicate  a  symmetric  boundary  condition  that  simulates  an 
infinite  plane. 

Above  and  below  the  unit  cell  are  regions  simulated  as  air.  These  regions  are  set 
up  using  a  Scattering  Boundary  Condition  with  Equation  (14)  as  the  boundary  conditions. 
On  the  upper  portion  of  the  model,  this  boundary  condition  will  act  a  sink  and  not  allow 
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any  reflected  waves  to  be  reflected  back  into  the  model.  All  waves  that  progress  all  the 
way  through  the  model  will  continue  to  the  Scattering  Boundary  Condition  below  the  unit 
cell  where  again  the  air  acts  as  a  sink  and  will  not  allow  the  waves  to  reflect  back  into  the 
model. 

n  X  (V  X  E)  -  jkn  x  (E  x  n)  =  -n  x  (Eg  x  (14) 

In  order  to  account  for  the  portion  of  the  incoming  wave  that  is  absorbed  by  the 
medium  and  microspheres,  a  volume  integration  of  resistive  heat  losses  is  performed  on 
the  material  of  concern.  The  incoming  THz  wave  is  set  at  1  W,  which  makes  for  a 
convenient  estimation  of  the  percentage  of  the  power  absorbed  by  the  material  of  interest. 
This  forms  the  basis  of  the  models  that  will  be  used  to  compare  to  the  experimentally 
derived  results. 

In  an  attempt  to  minimize  the  complexity  of  the  model  analyzed  within 
COMSOL,  a  unit  cell  model  was  the  basis  for  all  modeling  performed  in  support  of  this 
research.  The  basis  for  the  unit  cell  used  was  the  model  from  Alveset  ah,  discussed 
previously.  Various  configurations  of  unit  cells  were  built  to  see  which  most  closely 
replicated  the  experimentally  obtained  results. 

One  important  aspect  of  modeling  in  COMSOL  is  the  generation  of  the  mesh.  The 
mesh  size  should  be  such  that  the  elements  are  no  larger  than  one-fifth  of  the  incoming 
wavelength  in  size.  The  wavelengths  analyzed  in  COMSOL  were  set  up  such  that  5  THz 
to  18  THz  were  analyzed.  This  equates  to  wavelengths  between  16.7  pm  to  60  pm,  thus 
the  element  lengths  should  be  between  3.34  pm  to  12  pm.  As  the  wavelengths  get 
smaller,  the  mesh  will  become  very  dense.  This  will  add  many  degrees  of  freedom  to  the 
models.  Since  controlling  the  mesh  size  is  an  important  aspect  in  the  amount  of  time  that 
it  takes  for  each  finite  element  model  to  run,  a  balance  needs  to  be  made  as  to  a  good  size 
for  the  maximum  size  of  the  elements  used  in  the  models.  The  mesh  density  varied  for 
each  model  to  strike  a  balance  between  a  dense  enough  mesh  and  a  reasonable 
computation  time. 
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1. 


Medium  and  Substrate  Models 


In  an  effort  to  ensure  the  modeling  methodology  used  in  this  researeh  was  correct, 
simple  models  consisting  of  a  small  “unit  cell”  of  various  substrates  and  media  were  built 
and  analyzed.  This  allowed  for  reasonable  estimations  of  the  real  and  imaginary  parts  of 
the  refractive  indices  for  the  silicon  wafer  substrates  and  the  various  media  used.  A 
model  was  built  for  the  liquid  tape  on  a  silicon  wafer  sample. 

The  model  consisting  of  the  liquid  tape  only  sample  was  run  with  varying 
thicknesses  of  liquid  tape.  The  range  of  thickness  was  12  pm  to  19  pm.  The  refractive 
index,  both  real  and  imaginary  portions,  of  the  liquid  tape  was  determined  by  matching 
the  experimentally  obtained  liquid  tape  data  to  simulated  values  in  the  Thin  Film 
Analysis.  The  estimated  values  used  for  the  refractive  index  of  the  liquid  tape  were 
n=  2.84343  andk=0.2121. 

The  results  of  varying  the  liquid  tape  thickness  as  well  as  the  experimental  data 
are  plotted  in  Figure  40.  Each  micron  change  in  thickness  has  a  significant  effect  on  the 
position  of  the  peak  absorption  values.  The  actual  thickness  of  the  layer  of  liquid  tape  on 
the  silicon  wafer  was  measured  using  the  KLA  Tencor  D-120  Profilometer.  The  average 
value  obtained  equated  to  a  layer  13  pm  thick.  This  number  was  based  on  the  average  of 
three  measurements,  all  of  which  varied  significantly.  Based  on  the  results  from 
COMSOL,  it  would  appear  that  the  peaks,  seen  in  Figure  41,  most  closely  align  when  the 
thickness  of  the  liquid  tape  is  16  pm.  This  value  was  used  for  the  model  containing  one 
aluminum  sphere  in  liquid  tape,  to  be  discussed  later.  The  results  also  show  that  the 
values  for  both  the  real  and  imaginary  parts  of  the  refractive  index  are  not  100%  accurate 
and  require  further  refinement. 
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Figure  40  Absorption  Characteristics  of  a  Thin  Layer  of  Liquid  Tape  with 
Thicknesses  Varying  from  12  pm  to  18  pm 
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Figure  41  Absorption  Characteristics  of  a  16  pm  Layer  of  Liquid  Tape 

2,  A  Single  Microsphere  in  a  Thin  Medium 

The  simplest  model  built  was  a  single  microsphere  in  a  thin  medium.  This  model 
contained  one  aluminum  sphere,  20  pm  in  diameter,  on  16  pm  of  liquid  tape.  The  space 
above  the  sphere  is  air,  as  well  as  the  space  below.  The  modified  version,  based  on  the 
liquid  tape  model  mentioned  previously,  has  one  aluminum  sphere  partially  submerged  in 
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the  liquid  tape.  Both  models  are  in  Figure  42.  The  thiekness  of  the  liquid  tape  remained  at 
16  pm.  The  results  for  both  models  are  plotted  with  the  experimentally  obtained  data  in 
Figure  43.  Both  of  the  COMSOL  models  have  absorption  values  well  below  the 
experimentally  derived  values.  This  is  likely  due  to  inaeeurate  values  for  the  refractive 
indices.  Another  source  of  difference  between  the  experimental  and  the  results  from  the 
model  could  be  due  to  the  scattering  effects  possible  in  the  experimental  results.  The 
overall  shape  appears  to  be  close  to  the  experimentally  derived  values,  so  using  a  value  of 
16  pm  is  likely  close  to  the  actual  value  of  the  liquid  tape  thickness. 


Figure  42  COMSOL  Models  of  One  20  pm  Aluminum  Sphere  in  16  pm  of 

Liquid  Tape 
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Figure  43  Experimental  and  COMSOL  Results  for  Absorption  in  Liquid  Tape 

with  20  pm  Diameter  Aluminum  Spheres 


3,  Various  Numbers  of  Aluminum  Spheres  in  Photoresist  Solution 

A  series  of  COMSOL  models  were  built  using  the  SPR955  positive  photoresist  as 
a  medium  for  the  partieles.  The  models  started  off  with  minimal  eomplexity  and  were 
incrementally  changed  to  see  the  effects  of  varying  the  parameters.  The  parameters  varied 
were  the  number  of  aluminum  particles,  the  configuration  of  the  aluminum  particles,  the 
refractive  indices  of  the  photoresist  and  the  aluminum  particles,  and  the  conductivity  of 
the  aluminum  particles. 

One  property  varied  through  this  set  of  models  the  treatment  of  the  aluminum 
spheres.  Some  of  the  models  were  analyzed  using  the  conductivity  (a)  of  aluminum  to 
solve  the  finite  element  model.  The  value  used  in  this  case  was  a=  3.8  x  10^  S/m.  To  use 
the  refractive  index  values  of  aluminum,  further  calculations  were  required.  Experimentally 
derived  values  from  Ordal  et  al.  were  used  to  derive  a  2nd  order  polynomial  equation  for  n 
and  k  that  is  based  on  frequency  [20].  The  equations  used  were: 


=  2.13/'-60.24/  +  467.28 

(15) 

=  1.23/" -42.38/ +  490.40 

(16) 

where /is  the  frequency  in  THz.  The  values  and  graphs  used  to  determine  these  equations 
are  contained  in  the  Appendix. 
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Both  the  real  and  imaginary  portions  of  the  refractive  index  of  the  photoresist 
were  calculated  using  the  Thin  Film  Analysis  tools.  Due  to  the  frequency  dependence  of 
the  n  and  k  values,  the  experimental  data  was  analyzed  in  0.5  THz  increments  and  plotted 
to  find  equations  that  represent  the  frequency  dependence.  The  frequency  dependence 
was  minimal  for  the  real  portion  of  refractive  index,  so  a  linear  fit  was  used.  The 
imaginary  portion  had  more  variance  and  a  third  order  polynomial  was  used.  The  derived 
equations  are  Equations  (17)  and  (18): 

=0.0035/  + 1.5341  (17) 

=-0.00003/  +0.00085/  -0.00825/ +  0.05096  (18) 

where  /  is  the  frequency  in  THz.  The  graphs  used  to  derive  these  values  were  discussed 
previously  and  can  be  seen  in  Figure  37  and  Figure  38. 

Multiple  sets  of  experimental  data  were  looked  at  to  determine  the  refractive 
index  of  the  photoresist  and  were  used  for  various  models  in  COMSOF.  Equations  (17) 
and  (18)  provided  results  in  COMSOF  that  were  closer  to  the  experimental  results.  There 
were  differences  in  the  experimental  data  due  to  some  variance  in  the  thickness  of  the 
SPR  across  the  sample  and  due  to  the  amount  of  water  vapor  that  had  evaporated  from 
the  photoresist  as  the  sample  was  analyzed  on  multiple  days. 

The  first  two  models  built  consisted  of  two  aluminum  spheres  in  a  vertical  line 
and  then  four  aluminum  spheres  aligned  in  one  horizontal  plane.  These  models  used  a 
conductivity  of  3.8  x  lO’  S/m  for  aluminum  and  Equations  (17)  and  (18)  for  refractive 
index  of  the  photoresist.  The  COMSOF  values  obtained  for  absorption  did  not  reflect 
anything  close  to  the  experimentally  derived  values  but  showed  that  increasing  the 
number  of  spheres  from  two  to  nine  spheres,  as  seen  in  Figure  44  brings  the  model  closer 
to  the  experimental  results. 
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Figure  44  Absorption  Characteristics  of  Two,  Four,  Five,  Six  and  Nine  20  pm 

Spheres 


Due  to  the  results  seen  in  Figure  48,  the  results  between  the  six  and  nine  sphere 
models  were  fairly  similar.  In  order  to  strike  a  balance  between  computing  power  and 
accuracy,  the  six  sphere  model  was  selected  to  be  analyzed  in  depth  and  can  be  seen  in 
Figure  45.  The  model  was  built  so  that  there  was  20  pm  between  spheres  and  10  pm 
between  the  spheres  and  the  boundaries. 
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Figure  45  Basic  COMSOL  Model  for  Six  Aluminum  Spheres  in  Photoresist 

Solution 

The  variables  analyzed  for  the  six  sphere  model  were  a  set  n  and  k  versus  two 
different  frequency  dependent  refractive  index  values  for  the  photoresist,  and  a  series  of 
progressively  larger  aluminum  spheres  ranging  from  20  pm  to  40  pm.  The  spacing 
between  the  particles  remained  constant  as  the  particle  size  changed.  Changing  the  n  and 
k  for  the  SPR  only  had  the  effect  of  slightly  shifting  the  peaks  of  the  graph.  The  overall 
trend  of  the  graph  was  the  same  and  can  be  seen  in  Figure  46. 
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Figure  46  Absorption  Characteristics  of  Six  Spheres  in  Photoresist  Utilizing 

Various  Refractive  Indices 

Changing  the  sphere  sizes  had  a  significant  effect  on  this  model.  As  the  spheres 
got  larger,  the  peak  absorption  occurred  in  a  lower  frequency  range  seen  in  Figure  47. 
The  larger  spheres  absorb  more  of  the  lower  frequency  THZ  wave.  When  a  composite 
graph  of  the  three  sphere  sizes  is  constructed,  as  seen  in  Figure  48,  the  model  is  close  to 
the  actual  data  from  8  to  10.5  THz  and  a  very  close  match  from  10.5  to  18  THz.  A 
cursory  analysis  of  an  SEM  image,  seen  in  Figure  49,  shows  that  there  are  particles  that 
are  smaller  than  5  pm  and  they  are  other  particles  that  are  larger  than  35  pm.  Thus,  using 
a  variety  of  sizes  of  aluminum  spheres  is  a  reasonable  assumption  based  on  the  SEM 
images  taken  of  the  aluminum  powder. 
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Figure  47  COMSOL  Results  for  the  Absorption  of  Six  Aluminum  Spheres  in 
Photoresist,  with  Spheres  Ranging  from  20  pm  to  40  pm 


100 


— Expenmental 

/ 

^  / 

,  «  -  -  Composite  Absorption  of  20, 30  and  40  urn  A) 

Spheres 


70 


6S 

5  6  7  8  9  10  11  12  13  14  15  16  17  18 

THr 


Figure  48  Composite  Graph  of  Absorption  Characteristics  of  20,  30  and  40  pm 

Aluminum  Spheres 
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Figure  49  500x  SEM  Image  of  Aluminum  Powder  Showing  a  Range  of  Particle 

Sizes 


By  changing  the  sphere  size  from  20  pm  to  40  pm  and  maintaining  the  spacing 
between  the  spheres,  the  volume  fraction  increases  from  approximately  3.3%  aluminum 
to  5.7%  aluminum.  Since  the  overall  shape  of  the  graphs  is  fairly  consistent  with  the 
exception  of  shifting  the  peaks,  the  change  in  volume  fraction  may  not  have  had  a  large 
impact  on  the  overall  absorption  in  the  models. 

Using  the  COMSOL  models,  it  was  also  determined  that  most  of  the  heating 
occurs  in  the  upper  layers  of  the  spheres.  Figure  50  depicts  the  resistive  heating  losses  for 
the  six  sphere  model  at  13  THz  and  the  nine  sphere  model  at  13.8  THz.  The  frequencies 
were  selected  for  each  case  due  to  the  fact  that  the  maximum  absorption  occurs  at  the 
selected  frequency  for  each  model.  Both  models  show  that  most  of  the  THz  wave  is  in  the 
upper  layer  of  spheres  due  to  the  heating  that  occurs  there.  The  bottom  layers  in  both 
examples  show  minimal  heating.  Since  the  models  accurately  represented  the 
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experimental  results  above  10  THz,  adding  more  than  three  layers  of  spheres  is  not  likely 
going  to  add  any  more  accuracy  to  the  models. 


•jo” 


*  1.2.10"  (W/m*) 


2.S 


as 


ai 


» i.2s.io"(W/m^) 


» 1 12.10’  (W/m*) 


(a) 


(b) 


Figure  50  (a)  Resistive  bleating  Losses  (in  W/m  )  for  Six  Spheres  at  13  THz 

(b)  Resistive  Heating  Losses  (in  W/m  )for  Nine  Spheres  at  13.8  THz 
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V.  CONCLUSIONS 


The  work  in  this  thesis  foeused  on  developing  experimental  methods  to  determine 
the  absorption  in  the  THz  range  of  the  eleetromagnetie  speetrum  of  various  nano-  and 
miero-sized  partieles  in  a  variety  of  media.  The  results  were  then  used  to  build  and  verify 
finite  element  models  in  COMSOL  that  mirrored  the  experimentally  derived  results. 

The  experimental  results  showed  that  epoxy  resin  exhibits  large  amounts  of 
absorption  between  5  and  18  THz.  This  made  it  more  difficult  to  test  the  absorption 
characteristics  of  the  CNTs  and  the  nickel  particles  in  epoxy.  The  liquid  tape  absorbed 
much  less  of  the  incoming  THz  wave  between  5  and  18  THz,  but  it  proved  to  be  a 
difficult  medium  to  perform  accurate  experiment  on.  The  SPR9555  positive  photoresist 
solution  did  absorb  more  THz  than  expected,  but  the  experimental  results  proved  useful 
in  the  development  of  finite  element  models. 

A  significant  amount  of  the  finite  element  modeling  focused  on  the  aluminum 
spheres  in  the  SPR955  positive  photoresist.  The  research  showed  that  increasing  the 
number  of  aluminum  particles  in  a  unit  cell  to  six  spheres  was  an  accurate  representation 
of  the  experimental  results  above  10  THz.  The  inaccuracy  in  the  model  is  likely  due  to  a 
few  different  factors.  One  is  that  there  are  likely  particles  and  clumps  of  particles  that 
were  larger  than  40  pm.  Additionally,  extremely  accurate  values  for  the  refractive  index 
of  the  photoresist  were  difficult  to  obtain. 

Further  research  in  this  area  can  be  conducted  to  further  refine  the  COMSOL 
models  that  currently  exist.  As  previously  stated,  the  weakness  is  in  determining  accurate 
n  and  k  values  for  the  media  used.  The  Fabry-Perot  effects  seen  in  the  experimentally 
derived  data  of  all  of  the  samples  that  used  silicon  as  a  substrate  made  it  very  difficult  to 
use  the  tools  available  for  curve  fitting  due  to  the  periods  of  oscillations  seen  on  the 
experimental  results.  Additionally,  the  fact  that  many  of  the  media  used  displayed  a 
frequency-dependence  for  both  the  real  and  imaginary  portions  of  refractive  index  added 
more  uncertainty  to  the  values  of  n  and  k  that  were  used  in  the  COMSOL  models.  A 
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further  refinement  of  the  tools  used  to  determine  refraetive  index  eould  greatly  enhance 
the  product  of  the  COMSOL  models 

Additional  research,  using  the  refined  COMSOL  models,  can  focus  on  aspects 
other  than  the  absorption  of  THz,  such  as  the  Joule  heating  effects.  A  greater  focus  on  the 
thermal  effects  of  THz  radiation  on  various  composites  will  allow  for  better 
understanding  of  the  physical  changes  in  the  composite.  The  heating  of  CNTs  or  micro¬ 
particles  can  have  an  impact  on  the  strength,  hardness  and  fatigue  characteristics  of  a 
composite  material. 

Still  further  research  can  be  done  in  trying  to  build  models  involving  CNTs.  The 
CNTs  present  a  unique  challenging  in  finite  element  modeling.  When  CNTs  are  used  in  a 
composite,  there  are  a  number  of  characteristics,  such  as  the  length,  diameter  and 
chirality  that  vary  for  each  individual  CNT.  An  efficient  way  to  build  and  simulate  a 
model  capable  of  accurately  representing  realistic  results  is  a  logical  next  step  for  this 
research. 
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APPENDIX.  DETERMINATION  OF  REFRACTIVE  INDEX  FOR 

ALUMINUM 


Figure  5 1  Values  and  Equations  Used  for  the  Real  Portion  of  the  Refractive 

Index  of  Aluminum 
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Figure  52  Values  and  Equations  Used  for  the  Imaginary  Portion  of  the 

Refractive  Index  of  Aluminum 
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